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Abstract 

We analyze the tagged particle diffusion for kinetically constrained 
models for glassy systems. We present a method, focusing on the Kob- 
Andersen model as an example, which allows to prove lower and upper 
bounds for the self diffusion coefficient Ds- This method leads to the 
exact density dependence of Ds, at high density, for models with 
finite defects and to prove diffusivity, Ds > 0, at any finite density for 
highly cooperative models. A more general outcome is that under very 
general assumptions one can exclude that a dynamical transition, like 
the one predicted by the Mode-Coupling-Theory of glasses, takes place 
at a finite temperature/chemical potential for systems of interacting 
particle on a lattice. 

1 Introduction 

Many physical systems, in particular glass forming liquids, display a very 
slow dynamics at low temperature/high density [1]. The laboratory glass 
transition corresponds to the temperature/density at which the structural 
relaxation timescale becomes larger than the experimental one (e.g. one 
hour). At this point the glass-forming liquid falls out of equilibrium and 
becomes an amorphous rigid material called glass. Thus, the laboratory 
glass transition is nothing else than a dynamical crossover and not a true 
dynamical transition. However, a natural question is whether this dramatic 
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increasing of the structural relaxation time is due to an underlying dynamical 
transition that takes place at lower temperature/higher density (but that is 
unreachable experimentally). Indeed different analytical approaches, in par- 
ticular the Mode-Coupling-Theory (MCT) of the glass transition [2] predict 
a dynamical arrest at a finite temperature and chemical potential at which 
the structural relaxation timescale diverges and the diffusion coefficient of a 
tagged particle, called self-diffusion coefficient Ds, vanishes. The interpre- 
tation of this transition is based on the cage effect: particles are trapped 
in the cage formed by their neighbors (however see also [3]). For particles 
interacting with a smooth potential it is widely accepted that this MCT 
transition describes at most a dynamical crossover. Instead in the case of 
potentials with a hard core part, in particular hard sphere systems, there is 
no agreement. Experiments on colloids [4], that can be indeed modeled by 
hard sphere systems, are very often interpreted as if a real dynamical transi- 
tion took place. However, in these systems the microscopic timescale is much 
larger than for the other glass forming liquids (approximatively nine order of 
magnitude), thus it could be argued that one is just looking to the very first 
part of the increasing of the relaxation timescale, here MCT indeed applies 
but the point at which the dynamical MCT crossover takes place is shifted 
toward unreachable experimental timescales. 

What can one say theoretically about the existence of this type of tran- 
sition called sometimes dynamical arrest? First let us focus on the case of 
interacting particle systems on a lattice. If the potential between the particle 
is short-range, the Hamiltonian is not singular and the only constraint is the 
hard core one, i.e. maximum one particle per site, then it has been proved by 
Spohn [5,6] that on long time and length scales the tagged particle performs 
a simple Brownian motion with a self-diffusion coefficient that is positive 
at any finite temperature and any finite chemical potential (these models 
are called in the mathematical literature symmetric exclusion processes with 
speed change and spin exchange dynamics; in the following we'll refer simply 
to them as RLG, i.e. reversible lattice gases). Actually, one needs also the 
system to be ergodic but physically this is just due to the fact that if there 
is a very large correlation length and/or a very large correlation time than, 
in principle, the tagged particle will enter in the Brownian motion regime 
on larger length and time scales. So at a critical point the tagged particle 
could take an infinite time to enter in the Brownian motion regime. However 
the lower bound on the self diffusion coefficient proved in [5, 6] is valid re- 
gardless of the existence of a critical point. Thus, even if a phase transition 
would take place at a temperature and a chemical potential /ic one gets 
Ym^T^T,■^,^^,,Ds{T,^l) > 0. 

However, the lattice models that are known to display a whole phe- 
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nomenology analogous to glass forming liquids and for which the existence 
of a dynamical arrest at finite temperature/chemical potential has been sug- 
gested, are the ones in which further hard constraints besides hard core, are 
imposed. Among the most studied ones there are the kinetically constrained 
lattice gases (KCLG) for which the jump rates of particles are different from 
zero not only if the constraint of having maximum one particle per site is 
verified but also if some additional constraint is verified, hence the name 
kinetically constrained (see [7] for a recent review). We will comment on the 
extension of our work to statically constrained models, like the Lattice Glass 
Models introduced in [8], and more general cases in the conclusion. 

Note that the presence of additional hard constraints can change the 
physical mechanism behind tagged particle diffusion quite a lot. As explained 
in [7] the KCLG are characterized at high density by the existence of "defects" 
that are the analog of vacancies for the RLG. A rather good understanding 
has been reached for models in which defects, that can freely move in an 
otherwise completely filled lattice, consist in a finite number, independent 
of the particle density, of vacancies. Instead, the cases in which the motion 
of "defects" involve a number of vacancies that diverges when the particle 
density approaches one, were much less understood. Indeed in the literature 
one can find numerical simulations suggesting a dynamical transition at a 
density less than one [7]. Recently, in collaboration with D.S. Fisher [9], we 
have analyzed one of the most studied model of this type, the Kob- Andersen 
model [10] proving that no dynamical transition takes place whatsoever and 
unveiling what is the mechanism inducing the slow dynamics. Our results 
and techniques apply also to the other highly cooperative KCLG [11]. 

In this paper we show how one can, under very general assumptions, 
generalize the proof of diffusivity, Ds > of [5, 6] to KCLG (this has been 
already done for simple KCLG with finite size defects in [12]). Our aim is 
twofold. First, we want to show that under very general assumptions (see 
next sections and conclusion) a dynamical arrest cannot take place at finite 
temperature/chemical potential even if there are hard constraints other than 
the hard core one. Second, we want to present a rigorous technique and a 
method that allows to obtain upper and lower bounds on the self-diffusion 
coefficient for KCLG (or more general interacting particle systems). The 
usual techniques applied to KCLG in order to obtain predictions on Ds, as 
diagrammatic resummation or approximate closure of exact equations (see 
[7,13]), are completely out of control in the high density regime. Although 
some of them work well in the intermediate density regime compared to the 
results of numerical simulations, they fail in general at high density and, 
often, predict a spurious dynamical transition. Thus, in this context it is 
particularly important to have rigorous methods that allow one to obtain 
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solid analytical predictions. We will focus on the Kob- Andersen model as an 
example but our method can be applied to also other KCLG or more general 
interacting particle systems. 

Let us finally comment on the continuum case. We are not aware of any 
proof of diffusivity at low temperature/high density in the case of Hamilto- 
nian dynamics. Instead in the case of interacting Brownian particles, even 
with hard core (a good model for colloidal systems), it has been recently 
shown that the self diffusion coefficient is positive at any finite temper- 
ature/chemical potential under very general assumptions [14]. The proof 
works for hard spheres, Lennard- Jones potential, etc. . . and it is, as our proof, 
a generalization of the Spohn's proof [5,6] for RLG. 

Organization of the paper 
We think that the technique we will make use of is not well known in the 
community working on KCLG therefore we have written this article in a 
detailed and self-contained way. The expert reader may skip section |3] and 
quickly go through section |31 

We introduce in detail the KCLG in section |21 and we explain in detail the 
KA model in sectional In section |3] we recall some probabilistic techniques 
that have been used to prove bounds on the self diffusion coefficient for the 
RLG [6]. In sectional focusing on the KA model on a triangular lattice 
as an example, we explain how one can obtain the high density behavior 
of the self diffusion coefficient for KCLG with "defects" formed by a finite 
number of vacancies. This section is useful to introduce the notation and 
the method before facing the difficult case of highly cooperative KCLG. In 
section ini focusing on the KA model on a square lattice, we explain how one 
can obtain strictly positive lower bounds in the case of highly cooperative 
KCLG. Finally, in[7|we present a final discussion of our results. 

2 Kinetically constrained lattice gases 

In the last twenty years there has been a growing interest in kinetically 
constrained lattice gases. These were introduced as models for supercooled 
liquids close to the glass transition [15] and nowadays they are also studied 
as paradigm for general glassy systems [7]. KCLG are (apparently) similar 
to RLG. As discussed before there is however an important difference in the 
choice of the jump rates of particles that are different from zero not only if the 
constraint of having maximum one particle per site is verified but also if some 
additional constraint is verified, hence the name kinetically constrained. This 
choice of jump rates was originally devised in order to mimic the cage effect, 
that might be at the heart of the glassy behavior and the slow dynamics of 
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glass forming liquids. Indeed a molecule in a dense liquid is typically trapped 
in a cage created by surrounding particles (see [16] for a visual experimental 
example) and this takes place in the regime of temperature and density at 
which the dynamics slows down dramatically. 

More specifically, kinetically constrained lattice gases are stochastic lat- 
tice gases with hard core exclusion, i.e. systems of particles on a lattice A 
with the constraint that on each site there can be at most one particle. A 
configuration is therefore defined by giving for each site x G A the occupation 
number rj^ G {0, 1}, which represents an empty or occupied site respectively. 
The dynamics is given by a continuous time Markov process on the config- 
uration space = {0, 1}'^' which consists of a sequence of particle jumps. 
A particle at site x attempts to jump to a different site y with a fixed rate 
Cx,y{v)j which in general depends both on {x,y} and on the configuration 
T) over the entire lattice. The discretized time version of the process is the 
following. At time t choose at random a particle, let x be its position, and a 
site y. At time t + dt, the particle has jumped from x to y with probability 
Cxy{v{t))j while with probability 1 — Cxy{f]{t)) the configuration has remained 
unchanged. In other words, the probability measure at time t, fit, can be 
derived by the initial measure /xq as 

f^t{v)= Yl exp (Ct) fioiv') (1) 

»)'e{o,i}lA| 

where C, the generator of the Markov process, is the operator which acts on 
local functions / : I^a — ^ M as 

^fiv)= E ^^Mifm-m) (2) 

{x,y}cA 

where we defined 

{rjy a z = X 
V. if ^ = 2/ (3) 
Vz ii zj^x,y 

The simplest model is the simple symmetric exclusion process, SSEP, in 
which c^^y^^irj) = r]x{l — rjy) + rjy{l — -q^) for nearest neighbors {x,y}, 
Cx,y{v) = otherwise. Therefore, only nearest neighbor jumps are allowed 
and there are no further kinetical constraints besides hard core. The defi- 
nition kinetically constrained refers more properly to models in which jump 
rates impose additional requirements in order for the nearest neighbor move 
to be allowed. In other words, the rate Cx^yirf) can be zero for some choices 
of the configuration rj and the couple {x, ?/} even if 77^ = 1 77^, = 0, thus 
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preventing the jump of a particle from site x to final empty site y. From 
the above definition it is immediate to see that dynamics preserves the 
number of particles, i.e. the hyperplanes with fixed number N of particles 
f^A,7V XlxeA^^' ~ ^} invariant under dynamics. Moreover, 

in general the rates are chosen in order to satisfy detailed balance w.r.t. 
the uniform measure i^A,Ar on such hyperplanes. In other words condition 
Cx,y{v) = Cy^xij]^^) is satisfied for any choice of 77 and the couple {x,y}. This 
implies that the generator is reversible with respect to and therefore 
i^A,7v is stationary^. Note that i^a,7v is nothing else than canonical measure 
with zero Hamiltonian, i.e. with this choice of the rules there are no static 
interactions beyond hard core and an equilibrium transition cannot occur. 
However, a priori it is possible that a dynamical ergodic/non-ergodic tran- 
sition occurs for some choices of the rules. To our knowledge, the KCLG 
which have been considered so far do not display such transition (see [9,11] 
for the proof that an ergodic/non-ergodic transition does not occur in some 
highly cooperative KCLG). On the other hand, such transition has been 
proved to occur for a kinetically constrained spin model, namely North-East 
model [17,18]. We shall discuss in the following possible forms of such tran- 
sition rates. However, we emphasize since now that the degeneracy of the 
rates implies that z/A,Ar is not the unique invariant measure, i.e. the system 
is not ergodic on VL\,m and this will have several consequences on dynamics 
inducing a very different behavior with respect to RLG case. 

3 Definition of the Kob- Andersen model 

The KA model is a kinetically constrained lattice model with jump rates 



namely a particle can move only if the hard core constraint is verified, as for 
the SSEP, and only if both before and after the move it has no more than 
m neighboring particles. If A is an hypercubic d— dimensional lattice m will 
take values only from to 2(i — 1 (different values of m define different KA 
models). Note that for m = 2d — 1 the simple symmetric exclusion case 

^Let ^{g, h) = X^rjeo M('7)ff(^)^(^)- ^ is reversible with respect to /x if, for any functions 
/ and g, equality fi{g,£f) = fi{f,Cg) holds. By a direct calculation it is possible to 
check that detailed balance implies reversibility with respect to i^a.tv, therefore the choice 
g{rj) = 1 Vry, implies ^J.{Cf) = V/. This, together with |^ implies that z/A.Af is invariant 
under time evolution. 




(4) 
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with Hamiltonian equal to zero is recovered. For future purposes it is useful 
to reformulate the rule in term of motion of vacancies. Indeed, as can be 
easily verified, the above definition corresponds to vacancies moving only if 
the initial and final sites have at least s = z — m — 1 neighboring vacancies, 
with z = 2d the coordination number of the lattice. Therefore the model is 
completely defined by the choice of the couple d, m or equivalently d, s. 

Note that these rates satisfy detailed balance with respect to i^a^at, i.e. 
uniform measure on the hyperplanes with fixed number of particles. However 
there exist configurations that are blocked under the dynamics, is not the 
unique invariant measure on the hyperplane and the process is never ergodic. 
For example, in the case d = 2 s = 1 with periodic boundary condition, a 
configuration which has a double row of sites completely filled belongs to a 
different ergodic component with respect to any other configuration which 
does not contain such structure. Indeed, one can directly check that the 
particles belonging to the double row can never move. 

On the infinite lattice A = Z*^ the model satisfies detailed balance with 
respect to the Bernoulli product measure Hid p at any density p. It is im- 
mediate to check that for s > d the system is not ergodic at any density 
p > 0. Indeed in this case all d-dimensional hypercubes of any size which 
are completely occupied by particles are blocked forever. On the other hand, 
in [9, 11], we have proved that ioi s < d the model is ergodic at any density 
p e [0, 1] with probability one. As a by product of the ergodicity proof we 
have established the following property, which will be a key ingredient for the 
proof in sectional Let p be a positive number such that p < 1. For any fixed 
density p < 1 there exists a finite length S(p) such that, by sorting at random 
with probability p\^p a configurations on a hypercubic d-dimensional lattice 
A of linear size S(p), with probability greater than p this configuration is 
such that any particle exchange inside a finite box around the origin can be 
performed through a suitable path of allowed moves^. Furthermore p can be 
taken arbitrary close to one taking a suitable (large but finite) S(p). In the 
following we let a configuration with such property be a frameable conGg- 
uration and refer to [9,11] for the demonstration of above property and an 
explanation of the chosen name. 

4 Diffusion of the tagged particle 

Consider a kinetically constrained model on the infinite lattice A = Z*^ and 
start at time zero from the equilibrium distribution, so that the process 

^The length S(p) depends on the choice of d, s and the size of the box. we drop the 
dependence on these parameters for simphcity of notation. 
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will be stationary. Then single out one particle, the tracer, and analyze its 
motion. In the density regime where the process is ergodic one can repeat the 
arguments in [6,19] and show that under a diffusive rescaling the position of 
the tracer at time t, x(t), converges to a Brownian motion with self diffusion 
matrix Ds{p). More precisely, \\m^^Qex{e~'^t) = y/2Ds b{t), where b(t) is 
standard Brownian motion and the self diffusion matrix Ds{p) is given by 
the variational formula [5]'^: 

1 

{x,y}cA 

(5) 

where / is a unit vector in Z*^, (a-b) is Euclidean scalar product, r^yf] is 
the configuration obtained translating the configuration t] of y, fip^o is the 
Bernoulli measure at density p conditioned to the existence of a particle in 
the origin and the infimum is over all real-valued local functions /, i.e. func- 
tions which depend on a finite number of occupation variables. For spatially 
isotropic systems, as the one we will deal with, the self diffusion matrix is 
usually proportional to identity, the proportionality coefficient being the so- 
called self diffusion coefficient. For future purposes it is useful to define the 
two sums in ^ as -Di(p, /) and -D2(p, /) 



(l-DsipyT) = inf 



fir-yV'') - m?) (6) 
D2{pJ) = \ J2 f^pA^^Av)[fm-f{ri)?) (7) 

{x,y}cA 
x^O,y^O 

It is immediate to notice that -Di(p, /) > and D2{p, /) > at any density 
p, and therefore Ds{p) > 0. However, for the tagged particle process to be 

•^Note that by the identity lO one obtain that the matrix Ds is positive, therefore the 
square root is a well defined function. 
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diffusive, Ds should be strictly positive. For the RLG case in any dimension 
d > 2 the result -D^ > at any density p < 1 holds. Since the proof [6] 
uses as a key ingredient the fact that jumps from occupied to neighboring 
empty sites are always allowed, it cannot be trivially extended to kinetically 
constrained models. We emphasize that this is not only a technical difficulty. 
Indeed, due to the presence of kinetic constraints, the physical mechanism 
behind the diffusion of the tagged particle may be rather different as we shall 
show in the next sections. 

In the following section we will prove that Ds > at any p < 1 for 
KA model on a triangular lattice. For this choice of the lattice there exist 
"defects" , namely clusters of two vacancies that can freely move through an 
otherwise totally filled lattice and such that the tagged particle can move 
whenever such "defects" passes by. Therefore these defects play the same 
role of the vacancies for the RLG case and the above mentioned proof of 
diffusivity ( [6], see also [5])) can be generalized. In the following we sketch 
this result in some detail since it can be extended to other models with finite 
"defects" and it is useful to introduce the notation and the method before 
analyzing the more difficult case of KA model on hypercubic lattices. Indeed, 
in this case it is not possible to construct finite cluster of vacancies that can 
freely move into an otherwise totally filled lattice and additional work is 
required to prove diffusivity. As explained in [9], in this case the mechanism 
behind diffusion is based on the existence of "quasi-defects" , namely cluster 
of vacancies which can move into typical regions of the system and allow the 
motion of a tagged particle (via a particular path of elementary moves) when 
they pass by. The size of such defects is density dependent and diverges for 
p — > 1, which makes the motion highly cooperative. 

5 A model with finite size defects: KA s = 1 
on a triangular lattice 

Consider the KA model with s = 1 on a triangular lattice A represented in 
figure d More precisely, the set of sites in A is the union of sites in a square 
lattice Ai and in its dual A2, i.e. the lattice obtained by displacing Ai of 
(ei + e2)/2 with ei = (1,0), 62 = (0,1). Furthermore, two sites {x,y} G 
A are nearest neighbors if x — y = ±ei or x — y = ±(ei + 62) /2 or else 
X — y = ±(— Ci + 62)72. As already noticed, rates (jH) reformulated in terms 
of vacancies correspond to the rule that a vacancy can move only if the 
initial and final sites have at least s neighboring vacancies, where s = 1 
in the case we are considering. In the triangular lattice two neighboring 
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Figure 1: The triangular lattice. 



sites always share a common third neighbor. Therefore any of two nearest 
neighbor vacancies can move to the common third neighbor. In other words 
a set of two neighboring vacancies, which we call "a defect", can be freely 
moved into an otherwise totally filled lattice, as can be immediately checked 
(see figure HJ. 

5.1 Heuristic arguments and upper bound on the self 
diffusion coefficient 

There is a simple heuristic argument based on the independent motion of 
these defects that leads to the correct density dependence of the self-diffusion 
coefficient at high density. Call pd the density of defects and the timescale 
on which defects move. The self- diffusion coefficient Dg of a tagged particle 
is expected to be proportional to the inverse of the time Tp on which each 
particle moves of one step. On the timescale Td the number of particles that 
have jumped is of the order V pd where V is the total number of sites. Thus 
we find 

''-^Vpd oc Vp. (8) 

Td 

As a consequence, at density close to one, we get Ds oc pdj^d- Note that 
we have assumed that the size of the defects do not change, in particular do 
not diverge in the limit p — ^ 1, otherwise the reasoning has to be changed 
slightly. Since for the KA model on the triangular lattice we are focusing 
on the defects are formed by two vacancies we find, in the limit p — >■ 1, 
Prf oc (1 — p)^ and Td oc 0(1), hence, Dg oc (1 — p)^ (similarly in the SSEP 
case a defect is just a vacancy and one gets Ds oc (1 — p)). 

In the following we will show how, using the existence of these defects, it 
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Figure 2: Sequence of moves which allow the displacement of a couple of 
neighboring vacancies. Circles denotes empty sites, filled dots stand for oc- 
cupied sites 

is possible to generalize the proof of the RLG case [6] and show that indeed 
Ds oc (1 — p)^. In principle our procedure is generalizable to all cases with 
defects having a size which doesn't diverge in the limit p 1. To obtain 
Ds oc (1 — p)^ we shall prove that when p is close enough to one -D5 is 
bounded from above and below respectively by Ku{l — pY and Kl{1 — pY 
where Ku^l are two positive constants (of course Ku > Kl). 

The proof of the upper bound is very easy: it consists just in choosing 
an appropriate test function /(r/) and evaluating the term in the parenthesis 
in the variational formula (0). Indeed consider the test function /o(?7) = ^70, 
which for each configuration is equal to the occupation number in the origin. 
Plugging it into ^ we find a term which is proportional to the probability 
to have a vacancy on a fixed site and at least another vacancy close to it 
(this is nothing else that the probability to have a defect). Thus we obtain 
Ds <Ku{l- pf where Ku is a suitable positive constant. 

5.2 Lower bound on the self-difFusion coefficient 

Now let us focus on the lower bound which is more involved. Instead of 
dealing directly with the variational formula (0), we define a proper auxiliary 
model and proceed in two steps. First, we establish that D^"^ > 0; second, 
we prove that Ds > cDg^^ with c a positive constant. More precisely, we 
will introduce an auxiliary process and prove that the diffusion coefficient in 
direction ei is positive, i.e. eiD'^^^ei > and eiDsei > eiD^^^ei. In an 
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Figure 3: Sites inside the closed line correspond, form left to right, to sets 
Rl and Rl- 



analogous way one can then introduce an auxiliary process to show the same 
inequalities in direction (e2 + ei), which completes the proof of the positivity 
of the self diffusion matrix. 

5.2.1 Construction of the auxiliary process 

Let us introduce some notation. Consider the following subsets of A 



^1 
Rt 



[x + ei, ^ — I 

|(x-ei, x+ 1 



(9) 



namely i?^ are the two couples of neighboring sites of x represented in figure 

We next define r/^'^'^ as the configuration obtained from j] by exchanging 
the occupation numbers in R}^ with the corresponding ones in i?^, namely 



Vx—ei 


if 


Z 


= X + ei 


Vx+ei 


if 


Z 


= X — ei 




if 


z 


= x+^ 




if 


z 


= x+ 


Vz 


if 


z 





(10) 



Finally, we introduce the events 

^^={r/G(] : r/, = OVzei?i} A.. := ^VJ Al (11) 

which contain configurations having the two neighboring sites in i?^ empty, 
namely configuration with a defect in i?^. 
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Let us denote by x the position of the tagged particle and define the 
auxihary process as follows. At time t = the tagged particle is at the origin 
and there is a "defect" in a neighboring couple of sites, namely x(0) = and 

r]{0) e A- (i.e. r/ei(0) = r/e2/2+ei/2(0) = or r/_ej(0) = rj. e,/2+e2/2{0) = 0). 

Then the process evolves as follows 

(i) The tagged particles can jump from x to x + ei. The rate of the jump is 

one ifr] e Al (i.e. r]^+ei = and r]^_ei/2+e2/2 = 0); 

(ii) The tagged particles can jump from x to x — ei. The rate of the jump 
is one if 77 e (i.e. r]^_ei = and ?7a;-ei/2+e2/2 = 0); 

(iii) Configuration 77 is exchanged to r/^'^'^ (i.e. the occupation variables of 
sites x + ei, a; + 61/2 + 62/2 are exchanged with the occupation variable 
in X — 61 and x — 6i/2 + 62/2, respectively). The rate of the exchange 
is one; 

(iv) All other moves are not allowed. 

With this definition of the rules it is immediate to check that, since the 
configuration at time zero is such that 77(0) G Ax(o), at any subsequent time 
t > condition 77 (t) e Ax(t) will hold. Moreover, since the defect can be 
moved from R], to and viceversa (move (iii) above), the jump of the 
tagged particle to any of the two neighbors in direction 61 can always be 
performed via a sequence of at most two moves. Therefore the self diffusion 
coefficient in direction 61 is strictly greater than zero at any density p < 1, 
namely eiDg^^ei > 0. This can be rigorously proved in the same way as for 
RLG (see [5]). 



5.2.2 Proof of the inequality between Ds and D'g^'^ 

Let us now turn to the second step, namely establishing inequality Ds > 
cDg^^, with 6 > 0. Since move (iii) for the auxiliary process is not allowed 
for KA, some work is required to establish such inequality. The basic idea is 
to show that all the moves allowed for the auxiliary process can be performed 
through a proper finite sequence of elementary nearest neighbors jumps which 
are allowed for the original model. If this is true then it is natural to expect 
that the above inequality among the diffusion coefficients can be rigorously 
established. Let us outline the proof. Consider the second term, D2, of the 
variational formula (0). The above choice of the rates yields for the auxiliary 
process 
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As already mentioned any couple of neighboring vacancies can be freely 
moved through the lattice using elementary moves allowed by KA, see figure 
121 In particular, thanks to the fact that the t] which enters in (I12|l have a 
couple of vacancies either in the set Rl or Rl (or both), by using the basic 
moves in figure El it is possible to construct a path of elementary nearest 
neighbor exchanges which have unit rate for KA model and connect rj to 
j^i,2,x_ Moreover, such path is independent on the configuration outside R]. 
and R^. In other words, by recalling definition and letting the exchange 
operator T;j;,y be 

T.,yV = V''' (13) 

the following equivalence holds 

n 

V'''^=n^-.+i'-.^ (14) 

i=i 

where n is the length of the above defined path and rj —* rj^i+^'^i is the 
elementary exchange which constitute the j-th move of the path. Therefore, 
term /(r/^'^'^) — /(?]) in ()12j) can be rewritten by a telescopic sum as 

n 

1=1 

where we have defined rji, . . .rjn a.s rji = r], rji = rj^^i''^ . Then, by using (fTHjl 
and Cauchy-Scharwz inequality we obtain 

n 

Since the path r]i, . . .rjn has been chosen in order that the exchanges in the 
right hand side are all elementary exchanges allowed for KA process, i.e. 
with unit rate, we can rewrite the above inequality by introducing in the 
right hand side the corresponding jump rates 

n 

j=l 

Inserting the above inequality in ()12p. using the change of variables rji ^ rj 
and the invariance of equilibrium measure under exchange of occupation 
numbers yields, for any real-valued local function /, 
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Dru.p) = T-^Tj^<^AM{fiv'''n-fiv)Y>o 



{x,y}cA 
d{x,y)=l 



where we let Af be the maximal number of times a single exchange has been 
used in the path. In other words, let J\fij be the number of times operator 
T*-' appears in expression IT^ we let A/" := ma.x^ijjj\fij where the maximum is 
taken over all the couples of neighboring We emphasize that A/" and 

n are independent on the choice of the configuration, since the path is fixed 
once for all and does not depend on the value of the configuration outside 
Ax- iFTom the above inequality and recalling definition (jSJ it is immediate 
to conclude that 

^2(p,/)>^^^^/^r(p,/) (19) 

for any p and /. On the other hand, since the rates of the moves for the 
tagged particle in the auxiliary process are always smaller or equal to the 
correspondent rates for KA (see moves (i) and (iii)), the following inequality 
trivially holds among the first term of the variational expression (0) for 
for the two processes 

D,{pJ) > p,,o{AnDnpJ) (20) 
for any p and /. Therefore 

Dsip) > ^^^^^ D'^^^p) > cp,,o{An D^^^ip) (21) 

with c a strictly positive constant and pp^^A^) = (1— p'^— 4p^(l — p)— 4p^(l — 
p)^), which is also strictly positive at any density p < 1. In particular, in the 
high density limit, Ppfi{A^) oc (1 — p)^, hence, we finally get Ds > Kl{^ ^ 
with Kl a positive constant. 



6 Self diffusion coefficient for highly cooper- 
ative KCLG 

In the following we analyze the s = 1 KA model on a square lattice. This is 
a case in which it is impossible to identify defects with a density independent 
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size that can freely move inside the lattice. Thus, the key ingredient for the 
diffusivity proof discussed in the previous section does not hold. This is not 
a simple technical difficulty but it is deeply related to the fact that diffusion 
takes place in a different way in this case. As it has been found in [9] by 
analytical and numerical arguments the diffusion takes place thanks to the 
cooperative motion of a number of vacancies that diverges approaching unit 
density. This leads to an extremely rapid decreasing of the self-diffusion 
coefficient that, however, remains positive until unit density. 
In the following we want to show how one can prove that indeed Ds remains 
positive for p < 1. Our procedure is very general and can be applied to the 
other highly cooperative KA models as well as to many other interacting 
lattice particle systems. Note that in this case, at variance with what done 
in the previous section, we will not discuss the upper bound on Ds- Mainly 
because we do not expect that the lower bound we establish does give the 
right density dependence (see [9,11]). 

The strategy of the proof is similar to the one discussed in the previous 
section: we introduce a proper auxiliary process such that Ds > cD^^^ with 
c > and then we prove that D^"^ > 0. The key physical ingredient that 
allows us to find that such an auxiliary process exists is that ( [9,11] and 
section ^ all the particle exchanges inside a finite box around the origin 
can be performed with a very high probability p through a suitable path in 
configuration space that involves particles at most at distance from the 
origin. In this case the restricted configuration in the sublattice A= of size H 
around the origin is called frameable (see section IH}. Furthermore p can be 
chosen arbitrary close to one taking a suitably large H(p). If this property 
is verified then it's possible to find an auxiliary process that maps onto a 
random walk in a random environment and has Dg^^ > 0. 

Let us explain the idea in more detail. Consider a configuration on the 
infinite lattice 1? sorted at random with Bernoulli measure at density p and 
focus on a sub-lattice As of linear size H(p). As recalled in section (jH)), 
we know that the restriction of the configuration to A= is frameable with 
probability almost one. Therefore in the initial configuration the tagged 
particle is with very high probability inside a frameable region of size H. 
Moreover, if one divides the infinite lattice in sub- lattices of linear size S, 
there exists with unit probability a percolating cluster of sub-lattices such 
that the initial configuration restricted to each sub-lattice is frameable. 

Thus, roughly speaking, if we define an auxiliary process such that: (1) 
the tagged particle can move if it is inside a frameable square, (2) during the 
dynamics frameable sublattices remain frameable and the tagged particle 
remains always inside the percolating cluster then we can reconstruct any 
move of such process through a finite sequence of moves allowed by KA 
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Figure 4: Sites inside the closed line correspond, form left to right an up 
down, to sets Rt\ R-\ i?f , Rf, R-\ Rt\ K\ Qt\ Q".\ Qt\ Q^J 

(indeed any nearest neighbor move in a frameable configuration of linear size 
S can be performed through a sequence of elementary moves allowed by KA 
rules and the length of such path is at most of order S^)and prove inequality 
Ds > cDg^^. Furthermore since in the auxiliary process the particle moves 
on a percolating giant cluster then Dg^^ > 0. 

6.1 Construction of the auxiliary diffusive process for 
the KA s = 1 on a square lattice 

Let us introduce some notation. Consider the following subsets of 

-Ri"^^^ := {(a;i±ei,X2), (xi±ei,X2 + 62)} 

^i"^^^ ■= {(a:i±ei,a;2), (xi±ei,X2 - 62)} 

RI^^'^ '■= {(a:i,a;2±e2), (xi + ei,X2±e2)} 

R^^'^^ ■= {(a;i,X2±e2), (xi - ei,X2±e2)} 

Qi^'^ := {(xi,X2 + e2)Ui?i±i)} 

Q(±2) := {{x^,x,-e,)UR^^'^} (22) 

i?^ are the eight possible couples of neighboring sites {y, z} such that y, z x 
and |y — a;| = 1 or |z — x| = 1 (see figure |3)), in other words one among y 
and z is neighboring site to x and the couple does not contain x; Ql are the 
four possible choices of three sites that, together with x, form a two by two 
square (see figure E}. We next define r]^^ '-^^ , for i G {±1, . . . , ±4} as the 
configuration obtained from rj by exchanging the occupation numbers in R^^ 
with the corresponding ones in 

(^Rt\R^^\ ._ / ^-T2ei if ^ G Rt' /^o^ 

\ Vz if z^Rt'UR-' 
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Then, for i G {±1, ±2}, we introduce the events 

^L:={^efi : v. = 0\/zeQl} (24) 
i.e. configurations having all the sites of set Ql. empty and 

Bl:={ven : r^z = l yzERlve K] (25) 
for i G {±1, ±2, ±3, ±4}, where 



jri .= 1^ e . 3A G 5 : i?^ C A, 5(i?^, ^A) > 3, 7/|a G J'a} (26) 

and S is the set of squares in 1? of linear size at most H. Here 9A is the 
boundary of square sublattice A, ?7|a is the restriction of a configuration to 
A, is the set of frameable configuration in A and ^(A, S) is the minimum 
over the Euclidean distance of all the couples {x, ?/} with x E A and y E B. 
In other words, B^. is the set of configurations in which a pair of sites adjacent 
to X (region i?^) is filled and is internal to a frameable square of linear size 
at most S. Let //(O) be the configuration and x(0) the position of the tagged 
particle at time zero, we define f/(0) as 

^(0). = (l «f «■ (27) 

y rjz otherwise 

The dynamics of the auxiliary process is chosen as follows: 

(i) The tagged particle can move from x to x + ci. The jump has rate one 

if ?7 G A^^ and f/(0) G B'^^ or r/ G A^'^ and f/(0)i3^^, zero otherwise; 

(ii) The tagged particle can move from x to x — ei. The jump has rate one 
if 7] G and f/(0) G B~^ or 77 G A~'^ and 7^(0) G zero otherwise; 

(iii) The tagged particle can move from x to x + €2- The jump has rate one 
if ?7 G A^^ and f/(0) G or 77 G and f/(0) G zero otherwise; 

(iv) The tagged particle can move from x to x — €2- The jump has rate one 
if ?7 G A^"^ and f/(0) G or 77 G and f/(0) G zero otherwise; 

(v) Configuration 77 can be transformed in 77^^ , namely the exchange 
of occupation variables in R^^ and can be performed. The move 
has rate one if 77 G A^^ and 7/(0) G B^^ or 7/ G and 7/(0) G 
zero otherwise. 
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(vi) Configuration rj can be transformed in rj^"^ '^"^ , namely the exchange 
of occupation variables in i?^^ and can be performed. The move 
has rate one if 77 G A^'^ and f/(0) G B'"^ oi t] e A~'^ and r/(0) G 
zero otherwise. 

(vii) Configuration 77 can be transformed in rj^i^'^'^^ , namely the exchange 
of occupation variables in R'^^ and can be performed. The move 
has rate one if 77 G A'^^ and f/(0) G or 77 G A'^'^ and 7/(0) G 
zero otherwise; 

(viii) Configuration rj can be transformed in r]^^ , namely the exchange 
of occupation variables in i?^^ and can be performed. The move 
has rate one if 7/ G A~^ and 7/(0) G or 7/ G and 7/(0) G 
zero otherwise; 

(ix) All other moves are not allowed. 

In the following we will show that the above choice of the rates is a suitable 
choice to perform the proof of diffusivity, since the auxiliary process has a 
positive diffusion coefficient and any move can be reconstructed by a finite 
sequence of elementary moves allowed by KA. 

Consider an initial configuration such that the tagged particle is inside a 
frameable square of size S and such that all the sites in at least one of the sets 
^ito?' ^^x{o) (^^^ figure IH) are empty , where x(0) is the position of the tagged 
particle (i.e. 7/(0) G A], for some i). Then, both conditions will hold at any 
subsequent time. Indeed, moves (i)-(iv) are such that the tagged particle 
remains always inside the empty two by two square. On the other hand, 
moves (v)-(viii) are devised in order that the only vacancies that are moved 
during the process are those which belong at time zero to this two by two 
square, therefore sublattices of size S that are frameable at time zero remain 
frameable at later times^. The fact that moves for the auxiliary process 
occur always inside frameable regions of size at most S implies that any move 
can be performed through a finite sequence of elementary moves allowed by 
KA. Indeed, by the properties of frameable configurations, any move inside 
a configuration of size H can be performed by a sequence of order O(S^) 
moves with positive rate for KA dynamics [9,11]. By using path arguments 
analogous to those used in section Elfor the triangular case, it is then possible 
to establish inequality Ds > cD^^^, with c positive. Let us shortly recall how 
this argument works and emphasize an important difference occurring in this 

■*More precisely, sublattices that are frameable in 7^(0) are frameable also for ri(t) at 
any t. 
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case. For the triangular case we have defined an auxihary process such that 
any move of the latter can be performed by a finite path of at most n nearest 
neighbor moves allowed for the considered model. Such path does not depend 
on the choice of the configuration. Then, we have rewritten each term of 
the variational formula ^ of Z)|"^ as a telescopic sum on the exchanges 
along this path. Finally, by using Cauchy Schwartz inequality, the fact that 
each possible move is used at most twice in the path and by performing an 
exchange of variables, we concluded that D'^^ < cDs with c > 0. In this 
case we can proceed analogously, with the length n of the paths at most 
S^. However, the path now depends on the whole configuration. Indeed the 
sequence of allowed moves that one has to do in order to perform a pair 
exchange depends on the position of the vacancies, see [9,11]. This yields 
in the inequality among D'^^ and an overall factor M = 2^ besides the 
factor due to n = S^. Let us explain in some detail this statement. 
With a path argument analogous to the one done before, we rewrite each 
term in D'^^ corresponding to the exchange of particles in i?^, as a 
telescopic sum over allowed exchanges for KA, namely 

i=o 

where Cx,y are the jump rate of KA model, r/o = V^Vi^ ■ ■ ■ ^Vn = V^"" '^"^ is 
the path of allowed elementary moves which connects rj to 77^^ '^"^ and such 
that rji = Vj^^Ti'^^ for a couple of nearest neighbors {xj, To obtain above 

inequality we used Cauchy Schwartz inequality and the fact that n < S^. In 
order to obtain from left and right hand side the terms which appear in the 
variational formula (0) for the auxiliary process and for KA, respectively, 
we should average inequality ()28p over the Bernoulli measure conditioned to 
have a particle in zero. As we already emphasized, the sequence xq . . . x„ of 
sites in which the exchange is done, depends on the positions of vacancies in 
configuration rj. Therefore, if we do the change of variable rji ^ r} m (pH|l 
and use the invariance of measure under exchange of variables, many different 
terms on the left can give rise to the same term on the right. Actually, the 
crucial thing to know is the following. To each configuration 77 for which 
the exchange is allowed by the auxiliary process, associate the correspondent 
path ?7o, . . . ,?7„ in configuration space^. Then, for each elementary nearest 
neighbor exchange e, denote by Me the number of different configurations 77 

^Of course there could be different sequences to perform the same move. However, one 
can always give a prescription associating one of them for any choice of r/ and any give 
exchange. 
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that use such exchange and let M = maXe Me- Therefore M is the overall 
factor coming from possible overcounting of configuration when going from 
the mean of the left hand side of p8|l to the terms in the variational formula 
(jni) for Ds- Physically M takes into account the most severe dynamical 
bottleneck in phase space (see [9]). Moreover, since each path is composed of 
moves internal to the frameable region of size £ < S which contains the tagged 
particle, M is for sure less or equal to the total number of configurations inside 
a square of size S, namely inequality M < 2^ holds. Therefore, we finally 
obtain, closely enough to unit density, 

Ds>c{l-pf ^ or (29) 

where c is a positive constant. The term (1 — p^)fip,Ei^) comes from the 
condition that the configuration at time zero should have the tagged particle 
with three vacancies around and be inside a frameable square of size at most 
S and Dg^^ is the diffusion coefficient of the auxiliary process subject to this 
condition. This ends the first part of the proof. In the following section we 
shall show that indeed Dg^^ > 0. 

6.2 Lower bound for the self diffusion coefficient of the 
auxiliary process 

Let us now prove that Z)^"^ > 0, i.e. that diffusivity holds for the auxiliary 
process. In this case, the mechanism which guarantees diffusivity is different 
from the one discussed in section El Analogously to the KA on a triangular 
lattice (sectional), the auxiliary process we have just introduced is such that 
if the tagged particle is at time zero in a two by two square of vacancies (i.e. 
//(O) G ^^(0) for some i) and inside a larger frameable square of size at most 
S, both conditions will be always fulfilled at later times. However, now it is 
not true that the tagged particle can always be moved in a chosen direction 
Cj through a proper path. For example, if we want to move it in direction 
ei this is possible only if 77 G .4^^ or 77 G A'^'^. Otherwise, if 77 G 
or ?7 G A-^, the move is allowed only if before one makes the exchange 
rj — > ?7^^ ,Rx~i ^ jjRx ,Rx~2^ respectively. However these exchanges of 
rectangles (which are the analogous of exchange rj — s> rj^-'^^^^+'^i for RLG) 
are not always allowed. Indeed (see rules (v)-(viii)) they have positive rate 
only if in the initial configuration the rectangular regions R~^, respectively 
R~^, do not contain vacancies and are inside a frameable square of size at 
most H. Note that the rate of such exchanges (i.e. the rate of the exchange 
f] — > 77"^^ ,Rx~i g^j^^ _^ jjRx ,Rx-2 conditioned to the fact that rj G A^^ or 
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Tj G A^"^ respectively) does not depend on the configuration r/, but is fixed 
once for all by the choice of the initial configuration ?7(0). In other words 
the choice of the initial configuration fixes the good rectangles that can be 
exchanged. This observation will allow us to map the motion of the two 
by two square of three vacancies plus tagged particle to a random walk in 
a random environment corresponding to the cluster of good rectangles. We 
emphasize that this cluster does not change during dynamics, therefore the 
randomness of the environment is quenched. Note that the probability Pg for 
a given rectangle to be good is greater than (the probability that both the 
sites inside the rectangle are occupied) multiplied for the probability that it 
is inside a frameable region of size at most S which is almost one (this is 
the key physical ingredient). Therefore in the high density regime Pg is well 
above the threshold of conventional site percolation. This implies that with 
unit probability the initial configuration has a percolating cluster of good 
rectangles. By using that above the percolation threshold random walk on 
random environment has a positive diffusion coefficient [20] we will therefore 
obtain that the diffusion coefficient of the auxiliary process is strictly positive. 
In the following we will sketch the proof of the above argument in some detail. 

Let ?7(o,o) = ^7(0) be the initial configuration. Let us define the following 
sequence of configurations ?7(m,n) for (m, n) G 1? 



V{m.+l,n) 



V{m,n+1) 



x,x+ei 
''l{m,n) 



if r/( 



m,n) 



_+2 p-2 



'""^^ if Vin.,n) e A-', m e Bt^ 



<„:,nf if Vim,n) G , r/(0) G Bt' 



L <^,nf" if ^(-,n) e A~\ m e B^^ 



(30) 



(31) 
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x,x—ei 



if V{m,n) e A-\ r/(0) e B-^ 



V{m-l,n) < 



(32) 



r/f^ „f if r/(^,„) G r/(0) G 



(m,n— 1) 



x,x—e2 



I {m,n) 



if V{m,n) G A-^ r/(0) G 
if V{^,n) e At\ r/(0) G 



(33) 



where x = x{m, n) is the position of the tagged particle in configuration 
V(m,n) (we drop the dependence of x on m and n to get a more readable 
notation). Note that given ?7(m,n) and x{m,n) for a couple {m,n), using the 
above definitions one can reconstruct the whole sequence. 
Let us define a Markov process on with generator G acting on functions 
/ : (m, n) — > M as 



Gf{m,n) 



^ A- (m, n) {f{m + i,n)- f{m, n)) + 
i=±i 

^ A- (m, n) (/(m, n + i) - /(m, n)) 



(34) 



j=±i 



namely a two dimensional random walk with rates A^^ for the jump in direc- 
tion icj. We can now chose these rates in order that 77(m(t),n(t)) = vi't) for the 
auxiliary process. More precisely, we chose the rates in order that for any 
function f{r]) the expectation value over the probability fit evoluted with 
the generator of the auxiliary process coincides with the expectation over 
the measure on {■m{t),n{t)) generated by By considering the dynamics 
of the auxiliary process and definition (jS^), one can directly check that the 
choice of X^^i which satisfies above requirement is the following 
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A+i(m,n) = I^+i(^/{m,n))lB+inB+2(mO)) + 

^^+2(?7(„^,„))Ig+inB+2(?7(o,o)) + 

l^-i(^(m,n))%i(^(0,0)) + 
Iyt-2(^(m,n))lB+2('7(0,0)) 

(35) 

Aii(m,n) = I^-i(^/(m,n))lB;rini3-2(mo)) + 

I^-2('7(m,n))Ii5-lnB-2('7(0,0)) + 
I^+l(^{m,n))lB-i(^(0,0)) + 
1^+2 ('7{m,n)) 16^2(^7(0,0)) 

(36) 

I^2('7(m,n))lB+3(?7(0,0)) + 

I^^2(?7(m,„))Ig+4(?7(o,o)) 

(37) 

A^i(m,n) = I^2(r/(„,„))Ig-3ne-4(r/(o,o)) + 

2^+i('7(n^>ri))lB-3(?7(0,0)) 

Iyt-i(^(m,n))lB-*(mO)) (38) 

Let us define also 

A+i(m,n) = le+^nB+^iVm) 

A+i(m,n) = ls+3nB+4vm) 
A+i(m,n) = lB-^nB~4vm) 

(39) 

and 

Gf{m,n) := ^ A^ (m, n) (/(m + n) - /(m, n)) + 
i=±i 

^ Ai(m,n) (/(m,n + 2) - /(m,n)) (40) 
j=±i 
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^From definitions (jH^ and (jH^ it is immediate to check that A!|_^ > 
Because the variational formula (0) implies that the self diffusion coefficient 
is a monotonic increasing function of the jump rates we find that Ds for 
the random walk process with rates A^^ is certainly larger than the one 
for the process with rates A!|_^. By recalling definition and results on 
crossover length, we know that for sufficiently high density the probability 
w.r.t. Bernoulli measure of events 

n B-', 
Bf n Bt\ 
B-J n B-\ 

is almost one, and therefore greater then threshold probability for conven- 
tional percolation on the square lattice. Hence, by using the result in [20] 
which establishes a central limit theorem for random walk in random envi- 
ronment when bond probability is greater than percolation threshold we find 
that Ds for the process with rates A^^ (and therefore Ds for the process with 
rates A!|_]^) is strictly positive.^ Moreover, since when m goes to m -|- 2 (n to 
n-\-2) the first (second) coordinate of the tagged particle position increases 
at least of one unit, inequahty 

/it {xi{tf + X2{tf) > {m{tf + n{tf) > ct (41) 

holds, where /if is the evoluted of initial measure fip^ under the auxiliary 
process and c a positive constant. This allows us to conclude that Dg^^ > 
at any p < 1 which, together with inequality ()29|) . implies Ds > for KA 
model. 

7 Conclusion 

In this work we have presented a general procedure, focusing on the KA 
model as an example, that allows one to prove lower and upper bounds on 
interacting particle systems on a lattice in the case of vanishing rates. This 

^Note that the percolation problem we consider is a site percolation problem in which 
the site probability is correlated over a finite distance equal to 2S(p) whereas in [20] it 
was analyzed only the case with independent bond probability. However we expect that 
their result can be generalized to our case. Furthermore from the physical point of view 
there is no doubt that Ds will be positive in our case of correlated site disorder whenever 
a giant cluster exists (except at the critical point). 
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is a generalization of the Spohn's proof [5,6] for RLG. In particular, focusing 
on KA s = 1 model on a triangular lattice as an example, we show how 
to obtain the exact density dependence of Ds in the high density limit for 
KCLG with finite size defects (see also [12]). Whereas for highly cooperative 
KCLG, i.e. when the size of the defects diverges approaching unit density, 
our procedure allows to prove diffusivity, Ds > 0, at any density smaller 
than one as we have shown for the KA s = 1 model on a square lattice. 
As we have stressed previously our method is completely general and can 
be applied to general (short-range) interacting particle systems on a lattice, 
including systems with a non-trivial equilibrium measure as for example the 
statically constrained models introduced in [8]. In order to be successful, it 
just needs one crucial key physical property. One has to know that, for a 
random equilibrium configuration, all the particle exchanges inside a finite 
box around a fixed site, say the origin, can be performed with a very high 
probability p through a suitable path in configuration space (allowed by the 
kinetic rules) that involves particles at most at distance H from the origin. 
Furthermore one needs also that p can be chosen arbitrary close to one taking 
a suitably large S. 

This property is certainly valid if, given two equilibrium configurations taken 
at random inside a box of size L, one can show that there exists with prob- 
ability p' at least a path in configuration space, allowed by the kinetic rules, 
that connects the two configurations and that YmvL^^p' = 1 ^. Note that is 
just equivalent to say that the system at hand is irreducible in the thermo- 
dynamic limit. More correctly, in the thermodynamic limit one irreducible 
component VL covers all the configuration space, i.e. an equilibrium config- 
uration taken at random belongs to Vt with probability one. Therefore our 
results, combined with the ones obtained for Brownian interacting particle 
systems [14], strongly suggest that the only case in which a dynamical ar- 
rest (at which Ds vanishes at a finite temperature/chemical potential) might 
happen is only when a irreducible-reducible transition takes place and that 
one has to identify the two types of transition. 

To our knowledge, none of the short-range interacting particle systems which 
have been considered so far (on a lattice as well as in the continuum) has 
been proved to display such reducibility transition in dimension larger than 
one. 

Finally, we want to stress again that the non- vanishing of the self-diffusion 
coefficient Ds does not imply that the structural relaxation time scale Tq, 

^Tlie boundary condition for the box has to be chosen such that it is the worst possible 
in order to find the path connecting two configurations. For example for the KA one would 
choose boundary conditions that are equivalent to embed the box in a completely filled 
lattice. 
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cannot diverge. In particular, as discussed previously, at a second order 
phase transition the structural relaxation time scale diverges whereas the 
self-diffusion coefficient stays finite. Thus, the decoupling between Ds and 
I/tq is a necessary condition for the existence of an ideal glass transition 
taking place at finite temperature and chemical potential. In experiments 
on fragile liquids [21] a decoupling is indeed observed between Ds and I/tq 
(more correctly between Ds and the viscosity rf) but not as strong as one 
would have very close to a phase transition, i.e. DsTa oc Tq. 

It is pleasure to thank Daniel S. Fisher for all the discussions that we had 
on the subject of this work and, more generally, on kinetically constrained 
models for glasses. 
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